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We review the theoretical status of 6 ^ s hadronic penguin decays in the Standard 
Model and beyond. We summarize the main theoretical tools to compute Branching 
Ratios and CP asymmetries for 5 — > s penguin dominated nonleptonic decays, and 
discuss the theoretical uncertainties in the prediction of time-dependent CP asymme- 
tries in this processes. We consider general aspects of 5 ^ s transitions beyond the 
Standard Model. Then we present detailed predictions in supersymmetric models 
with new sources of flavor and CP violation. 



I. INTRODUCTION 



New Physics (NP) can be searched for in two ways: either by raising the available energy 
at colliders to produce new particles and reveal them directly, or by increasing the experi- 
mental precision on certain processes involving Standard Model (SM) particles as external 
states. The latter option, indirect search for NP, should be pursued using processes that 
are forbidden, very rare or precisely calculable in the SM. In this respect. Flavor Changing 
Neutral Current (FCNC) and CP-violating processes are among the most powerful probes 
of NP, since in the SM they cannot arise at the tree-level and even at the loop level they 
are strongly suppressed by the GIM mechanism. Furthermore, in the quark sector they 
are all calculable in terms of the CKM matrix, and in particular of the parameters p and 
f] in the generalized Wolfenstein parametrization [1]. Unfortunately, in many cases a deep 
understanding of hadronic dynamics is required in order to be able to extract the relevant 
short- distance information from measured processes. Lattice QCD allows us to compute the 
necessary hadronic parameters in many processes, for example in AF = 2 amplitudes. In- 
deed, the Unitarity Triangle Analysis (UTA) with Lattice QCD input is extremely successful 
in determining p and f] and in constraining NP contributions to AF = 2 amplitudes [2-6]. 

Once the CKM matrix is precisely determined by means of the UTA (either within the 
SM or allowing for generic NP in AF = 2 processes), it is possible to search for NP con- 



2 



tributions to AF — 1 transitions. FCNC and CP-violating hadronic decays are indeed the 
most sensitive probes of NP contributions to penguin operators. In particular, penguin- 
dominated nonleptonic B decays can reveal the presence of NP in decay amplitudes [7-9] . 
The dominance of penguin operators is realized in 6 — sqq transitions. 

Thanks to the efforts of the BaBar and Belle collaborations, i?-factories have been able 
to measure CP violation in several 6 — > s penguin-dominated channels with an impressive 
accuracy [10-20]. To fully exploit this rich experimental information to test the SM and 
look for NP, we need to determine the SM predictions for each channel. As we shall see 
in the following, computing the uncertainty in the SM predictions is an extremely delicate 
task. Only in very few cases it is possible to control this uncertainty using only experimental 
data; in general, one has to use some dynamical information, either from flavor symmetries 
or from factorization. Computing CP violation in 6 — > s penguins beyond the SM is even 
harder: additional operators arise, and in many cases the dominant contribution is expected 
to come from new operators or from operators that are subdominant in the SM. In the 
near future, say before the start of the LHC, we can aim at establishing possible hints of 
NP in 6 ^ s penguins. With the advent of the LHC, two scenarios are possible. If new 
particles are revealed, b ^ s penguin decays will help us identify the flavor structure of the 
underlying NP model. If no new particles are seen, 6 — > s penguins can either indirectly 
reveal the presence of NP, if the present hints are confirmed, or allow us to push further the 
lower bound on the scale of NP. In all cases, experimental and theoretical progress in 6 — > s 
hadronic penguins is crucial for our understanding of fiavor physics beyond the SM. 

This review is organized as follows. In Sec. II we quickly review the basic formahsm for 
b ^ s nonleptonic decays, and the different approaches to the calculation of decay amplitudes 
present in the literature. In Sec. Ill, we present the predictions for Branching Ratios (BR's) 
and CP violation within the SM following the various approaches, and compare them with 
the experimental data. In Sec. IV, we discuss the possible sources of NP contributions to 
b ^ s penguins and how these NP contributions are constrained by experimental data on 
other b ^ s transitions. In Sec. V, we concentrate on SUSY extensions of the SM, discuss 
the present constraints and present detailed predictions for CP violation in 6 ^ s penguins. 
In Sec. VI we briefly discuss 6 — > s penguins in the context of non-SUSY extensions of the 
SM. Finally, in Sec. VII we summarize the present status and discuss future prospects. 
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II. BASIC FORMALISM 



A. 



Generalities 



The basic theoretical framework for non-leptonic B decays is based on the Operator 
Product Expansion (OPE) and renormahzation group methods which allow to write the 
amphtude for a decay of a given meson B—B^, Bg, -B+ into a final state F generally as 
follows: 



Here Tieff is the effective weak Hamiltonian, with Qi denoting the relevant local operators 
which govern the decays in question within the SM, and Qi denoting the ones possibly 
arising beyond the SM. The CKM factors and the Wilson coefficients Ci{fi) describe 

the strength with which a given operator enters the Hamiltonian; for NP contributions, we 
denote with Cf^{ii) and Cf^{ii) the Wilson coefficients arising within a given NP model, 
which can in general be complex. In a more intuitive language, the operators Qi{ij) can 
be regarded as effective vertices and the coefficients Ci(//) as the corresponding effective 
couplings. The latter can be calculated in renormalization-group improved perturbation 
theory and are known including Next-to-Leading order (NLO) QCD corrections within the 
SM and in a few SUSY models [21-23]. The scale /j, separates the contributions to A{B F) 
into short-distance contributions with energy scales higher than /i contained in Cj(/x) and 
long-distance contributions with energy scales lower than /i contained in the hadronic matrix 
elements {Qi{lJ'))- The scale /i is usually chosen to be 0(mb) but is otherwise arbitrary. 
The effective weak Hamiltonian for non-leptonic 6 — > s decays within the SM is given by: 




A^NP 



+ E^r(/^)(^l4(/^)|i?). 



(1) 



i=l 
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with 



Q3,5 — 

Qn = 



{hLYu\){u\lnSL), 
Y.q(hLl^SL){qL,Rl^iqL,R) , 
-'^HqipLqRjiqRSL) , 

-3 Eg eq(bLqR){qRSL) , 



Q4 = Eq(bLl^qL)(qLl^cSL) , 



Q7,9 = i llq{hLl^SL)eq{qLl^qL) , 

Qio = i Eqey{bLi''qL){qLit.sL) , 



(3) 



where g^^j? = (1 ^ ^^)/2q, — {«, c} and denotes the quark electric charge (e^ = 2/3, 
= — 1/3, etc.). The sum over the quarks q runs over the active flavors at the scale /x. 
Qi and Q2 are the so-called current-current operators, Qz-q the QCD-penguin operators, 
Q7-10 the electroweak penguin operators and (5ii,i2 the (chromo)-magnetic penguin oper- 
ators. C.t(/i) arc the Wilson coefficients evaluated at = ©(m^). They depend generally 
on the renormalization scheme for the operators. The scale and scheme dependence of the 
coefficients is canceled by the analogous dependence in the matrix elements. It is therefore 
convenient to identify the basic renormalization group invariant parameters (RGI's) and to 
express the decay amphtudes in terms of RGI's. This exercise was performed in ref. [24], 
where the RGI's were identified and the decay amplitudes for several two-body nonleptonic 
B decays were written down. For our purpose, we just need to recall a few basic facts about 
the classification of RGI's. First of all, we have six non-penguin parameters, containing only 
non-penguin contractions of the current-current operators Qi,2'- emission parameters £^1,2, 
annihilation parameters Ai^2 and Zweig-suppressed emission- annihilation parameters EAi^2- 
Then, we have four parameters containing only penguin contractions of the current-current 
operators (5i,2 in the GIM-suppressed combination (55,2 ~ Qi,2- P^^^ and Zweig suppressed 
P^^^. Finally, we have four parameters containing penguin contractions of current-current 
operators (5i,2 (the so-called charming penguins [25]) and all possible contractions of penguin 
operators (53-12 ^ Pi,2 and the Zweig-suppressed Ps ^. 

Let us now discuss some important aspects of 6 ^ s penguin nonleptonic decays. First 
of all, we define as pure penguin channels the ones that are generated only by Pj and pP^^ 
parameters. Pure penguin 6 — > s decays can be written schematically as: 



A{B^F)^-V:,VusY.Pl 



iGIM 



(4) 



Neglecting doubly Cabibbo suppressed terms, the decay amplitude has vanishing weak phase. 
Therefore, there is no direct CP violation and the coefficient Sp of the sin Ami term in the 
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time-dependent CP asymmetry (for F a CP eigenstate with eigenvalue 77^?) measures the 
phase of the mixing amphtude: Sp — r^pImA^ = —r]Fsm2(f)M, where Xp = — e~'^^'f''^ , 
A = A{B ^ F), A = A{B F) and 4>m = (3 (-A) for Bd (B.,) mixing. Comparing the 
measured Sp to the one obtained from b ccs transitions such as B^^s) J/'^Ks{(f)) can 
reveal the presence of NP in the b ^ s penguin amplitude. However, to perform a precise 
test of the SM we need to take into account also the doubly Cabibbo suppressed terms in 
Eq. (4). The second term then acquires a small and calculable weak phase, leading to a small 
and calculable AS" = —tifSf — sin 20m- Furthermore, we must consider the contribution 
from the first term, i.e. the contribution of GIM penguins. An estimate of the latter requires 
some knowledge of penguin-type hadronic matrix elements, which can be obtained either 
from theory or from experimental data. Let us define this as the "GIM-penguin problem" : 
we shall come back to it in the next Section after introducing the necessary theoretical 
ingredients. 

Besides pure penguins, we have 6 — > s transitions in which emission, annihilation or 
emission-annihilation parameters give a contribution to the decay amplitude. Let us call 
these channels penguin-dominated. Then we can write schematically the decay amplitude 
as: 

A{B ^ F) = -V:,V^s E {Ti + i^^^^) - V*Vts E P. , (5) 

where Tj = {Ei,Ai,EAi}. Also in this case, neglecting doubly Cabibbo suppressed terms 
the decay amplitude has vanishing weak phase, so that AS = at this order. However, we 
expect Ti > Pj so that the double Cabibbo suppression can be overcome by the enhancement 
in the matrix element, leading to a sizable A^*. Once again, the evaluation of AS" requires 
some knowledge of hadronic dynamics. Let us define this as the "tree problem" and return 
to it in the next Section. 

B. Evaluation of hadronic matrix elements 

The last decade has witnessed remarkable progress in the theory of nonleptonic B decays. 
Bjorken's color transparency argument has been put on firm grounds, and there is now a wide 
consensus that many B two-body decay amplitudes factorize in the limit mb — > 00 and are 
therefore computable in this limit in terms of few fundamental nonperturbative quantities. 
Three different approaches to factorization in B decays have been put forward: the so-called 
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QCD factorization [26-30], perturbative QCD (PQCD) [31-33] and Soft-CoUinear Effective 
Theory (SCET) [34-38]. A detailed discussion of these approaches goes beyond the scope 
of this review; for our purpose, it suffices to quickly describe a few aspects that are relevant 
for the study of 6 ^ s penguin nonleptonic decays. 

Unfortunately, as suggested in ref. [25] and later confirmed in refs. [39-58], it turns out 
that subleading corrections to the infinite mass limit, being doubly Cabibbo-enhanced in 
6 — > s penguins, are very important (if not dominant) in these channels, so that they 
reintroduce the strong model dependence that we hoped to eliminate using factorization 
theorems. While different approaches to factorization point to different sources of large 
corrections, no approach is able to compute from first principles all the ingredients needed 
to test the SM in 6 ^ s penguins. Therefore, it is important to pursue, in addition to 
factorization studies, alternative data-driven approaches that can in some cases lead to 
model-independent predictions for CP violation in & — > s penguins. 

Let us now quickly review the main tools that are available for the study of 6 — > s 
penguins. 



The first step towards a factorization theorem was given by Bjorken's color transparency 
argument [59]. Let us consider a decay of the B meson in two light pseudoscalars, where 
two light quarks are emitted from the weak interaction vertex as a fast-traveling small-size 
color-singlet object. In the heavy-quark limit, soft gluons cannot resolve this color dipole 
and therefore soft gluon exchange between the two light mesons decouples at lowest order 
in A/rrib (here and in the following A denotes a typical hadronic scale of order Aqcd)- 

Assuming that in B decays to two light pseudoscalars perturbative Sudakov suppression is 
not sufficient to guarantee the dominance of hard spectator interactions, QCD Factorization 
(QCDF) states that all soft spectator interactions can be absorbed in the heavy-to-hght form 
factor [26]. Considering for example S — > tttt decays, the following factorization formula 
holds at lowest order in A/nib'- 



1. QCD factorization 



Jo 



(6) 
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where f^^'^{q^) is a, B ^ n form factor, and (j)-^ {4>b) are leading- twist light cone distribution 
amphtudes of the pion (B meson). T/'^^ denote the hard scattering amplitudes. Notice that 
starts at zeroth order in as and at higher order contains hard gluon exchange not involving 

the spectator, while T^^ contains the hard interactions of the spectator and starts at order 
as- 

The scheme and scale dependence of the scattering kernels j/'^^ matches the one of 
Wilson coefficients, and the final result is consistently scale and scheme independent. 

Final state interaction phases appear in this formalism as imaginary parts of the scattering 
kernels (at lowest order in A/rrib). These phases appear in the computation of penguin 

contractions and of hard ghion exchange between the two pions. This means that in the 
heavy quark limit final state interactions can be determined perturbatively. 

A few remarks are important for the discussion of CP violation in b ^ s penguins: 

• Penguin contractions (including charming and GIM penguins) are found to be factor- 
izable, at least at one loop. 

• Subleading terms in the A / rrih expansion are in general non-factorizable, so that they 
cannot be computed from first principles. They are important for phenomenology 
whenever they are chirally or Cabibbo enhanced. In particular, they cannot be ne- 
glected in b ^ s penguin modes. This introduces a strong model dependence in the 
evaluation of 6 — > s penguin BR's and CP asymmetries. 

• Power suppressed terms can invalidate the perturbative calculation of strong phases 
performed in the infinite mass limit. Indeed, in this case subleading terms in the 
A/mb expansion can dominate over the loop-suppressed perturbative phases arising at 
leading order in A/m^. 

2. PQCD 

The basic idea underlying PQCD calculations is that the dominant process is hard gluon 
exchange involving the spectator quark. PQCD adopts the three-scale factorization theorem 
[60] based on the perturbative QCD formahsm by Brodsky and Lepage [61], with the inclu- 
sion of the transverse momentum carried by partons inside the meson. The three different 
scales are the electroweak scale M^, the scale of hard gluon exchange t ~ 0(-\/Am^), and 
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the factorization scale 1/6, where b is the conjugate variable of parton transverse momenta. 
The nonperturbative physics at scales below 1/6 is encoded in process-independent meson 
wave functions. The inclusion of transverse momentum leads to a Sudakov form factor which 
suppresses the long distance contributions in the large b region, and vanishes as 6 > 1/A. 
This suppression renders the transverse momentum flowing into the hard amplitudes of order 
Arrib. The off-shellness of internal particles then remains of 0{Amb) even in the end-point 
region, and the singularities are removed. 
Notice that: 

• Contrary to QCD factorization, in PQCD all contributions are assumed to be cal- 
culable in perturbation theory due to the Sudakov suppression. This item remains 
controversial (see refs. [62] and [63]). 

• The dominant strong phases in this approach come from factorized annihilation dia- 
grams. 

• Also in this case, there is no control over subleading contributions in the A/m^ expan- 
sion. 

3. SCET 

Soft-coUinear effective theory is a powerful tool to study factorization in multi-scale prob- 
lems. The idea is to perform a two-step matching procedure at the hard {0{mi,)) and hard- 
coUinear {0{\/rnJi)) scales. The final expression is given in terms of perturbative hard 
kernels, light-cone wave functions and jet functions. For phenomenology, it is convenient to 
fit directly the nonperturbative parameters on data using the following expression for the 
decay amplitude, valid at leading order in as [56-58]: 

A{B ^ M1M2) oc /m.Cj""^ /' du4>MAu)Tu{u) + fM.C^'T,^ + 1^2 + A^^^^ , 

JO 

where T's are perturbative hard kernels, C's are nonperturbative parameters and Acc denotes 
the "charming penguin" contribution. 
We notice that: 

• Charming penguins are not factorized in the infinite mass limit in this approach, 
contrary to what obtained in QCD factorization; 
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• Phenomenological analyses are carried out at leading order in ag and at leading power 
in A/mb, 

• No control is possible on power corrections to factorization. 

4- SU{3) flavor symmetry 

An alternative approach that has been pursued extensively in the literature is to use 
SU{3) flavor symmetry to extract hadronic matrix elements from experimental data and 
then use them to predict S'[/(3)-related channels [64-87]. In principle, in this way it is 
possible to eliminate all the uncertainties connected to factorization and the inflnite mass 
limit. On the other hand, SU (3)-breaking must be evaluated to obtain reliable predictions. 

A few comments are in order: 

• In some fortunate cases, such as the contribution of clcctroweak penguins Qg.io to 
B Kn decays, SU (3) predicts some matrix elements to vanish, so that they can be 
assumed to be suppressed even in the presence of 5'C/(3) breaking [88-90]. 

• Explicit nonperturbativc calculations of two-body nonleptonic B decays indicate that 
SU (3)-breaking corrections to B decay amplitudes can be up to 80%, thus invalidating 
SU{?>) analyses of these processes [91]. 

• To take partially into account the effects of S'f/(3) breaking, several authors assume 
that symmetry breaking follows the pattern of factorized matrix elements. While this 
is certainly an interesting idea, its validity for 6 — > s penguins is questionable, given 
the importance of nonfactorizable contributions in these channels. 

5. General parameterizations 

The idea developed in Rcfs. [42, 92] is to write down the RGI parameters as the sum 
of their expression in the infinite mass limit, for example using QCD factorization, plus an 
arbitrary contribution corresponding to subleading terms in the power expansion. These 
additional contributions are then determined by a fit to the experimental data. In 6 — > s 
penguins, the dominant power-suppressed correction is given by charming penguins, and the 
corresponding parameter can be determined with high precision from data and is found to 
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be compatible with a A/ nib correction to factorization [42]. However, non-dominant correc- 
tions, for example GIM penguin parameters in 6 ^ s decays, can be extracted from data 
only in a few cases (for example in B ^ Kit decays) [92]. However, predictions for AS 
depend crucially on these corrections, so that one needs external input to constrain them. 
One interesting avenue is to extract the support of GIM penguins from SU (3)-related chan- 
nels {b ^ d penguins) in which they are not Gabibbo-suppressed, and to use this support, 
including a possible SU{3) breaking of 100%, in the fit of 6 — > s penguin decays. Alter- 
natively, one can omit the calculation in factorization and fit directly the RGI parameters 
from the experimental data, instead of fitting the power-suppressed corrections [93, 94]. 
We remark that: 

• Gompared to factorization approaches, general parameterizations have less predictive 
power but are more general and thus best suited to search for NP in a conservative 
way. 

• This method has the advantage that for several channels, to be discussed below, the 
predicted AS decreases with the experimental uncertainty in BR^s and GP asymme- 
tries of 6 — > s and SU (3)-related b ^ d penguins. 

We conclude this Section by remarking once again that neither the "GIM-penguin prob- 
lem" nor the "tree problem" can be solved from first principles and we must cope with 
model-dependent estimates. It then becomes very important to be able to study a variety 
of channels in several different approaches. In this way, we can hope to be able to make 
solid predictions and to test them with high accuracy. In the following, we quickly review 
the present theoretical and experimental results, keeping in mind the goal of testing the SM 
and looking for NP. 

III. BR'S AND CP ASYMMETRIES WITHIN THE SM 

The aim of this Section is to collect pre- and post-dictions for BR^s and GP asymmetries 
of 6 — > s penguin decays obtained in the approaches briefly discussed in the previous Section. 
The main focus will be on AS, but BR^s and rate GP asymmetries will play a key role in 
assessing the reliability and the theoretical uncertainty of the different approaches. 
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A. BR^s and rate CP asymmetries 

In Tables I-III we report some of the results obtained in the literature for B decay BR^s 
and CP asymmetries. For QCD Factorization (QCDF) results, the first error corresponds to 
variations of CKM parameters, the second to variations of the renormalization scale, quark 
masses, decay constants (except for transverse ones), form factors, and the rj — rj' mixing 
angle. The third error corresponds to the uncertainty due to the Gcgcnbauer moments in the 
expansion of the light-cone distribution amplitudes, and also includes the scale-dependent 
transverse decay constants for vector mesons. Finally, the last error corresponds to an 
estimate of the effect of the dominant power corrections. For PQCD results from refs. [54, 
55], the error only includes the variation of Gegenbauer moments, of |T4b| and of the CKM 
phase. For PQCD results from ref. [95], the errors correspond to input hadronic parameters, 
to scale dependence, and to CKM parameters respectively. For SCET results, the analysis is 
carried out at leading order in and A/ nib assuming exact SU (3). The errors are estimates 
of SU{3) breaking, of A/rrib corrections and of the uncertainty due to SCET parameters 
respectively. SCET I and SCET II denote two possible solutions for SCET parameters in 
the fit [58]. For General Parametrization (CP) results, the errors include the uncertainty 
on CKM parameters, on form factors, quark masses and meson decay constants, and a 
variation of A/m^ corrections up to 50% of the leading power emission amplitude. The 
values in boldface correspond to predictions {i.e. the experimental value has not been used 
in the fit). 

First of all, we notice that all approaches are able to reproduce the experimental BR^s 
of S — > PP penguins, although QCDF tends to predict lower BR^s for B — > Prj', albeit 
with large uncertainties. Concerning BR^s oi B ^ PV penguins, QCDF is always on the 
low side and reproduces experimental BR^s only when the upper range of the error due to 
power corrections is considered. PQCD shows similar features for K* and p modes, while it 
predicts much larger values for BR^s of B ^ Kuj decays. 

The situation for rate CP asymmetries is a bit different. Both QCDF and SCET predict 
AcY>{B^ — > 7r°K°) ~ —AcY>{B^ — > Tr+X") while experimentally the two asymmetries have 
the same sign. PQCD reproduces the experimental values, although it predicts Acp(-B° — > 
TT^K^) on the low side of the experimental value. It is interesting to notice that the GP 
approach is able to predict the correct value and sign of Acp{B^ —>■ tt^K^) in spite of the 
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TABLE I: Results for CP-averaged BR^s (in units of 10^^) and CP asymmetries (in %) in several 
approaches for B — > PP decays. Experimental averages from the Heavy Flavor Averaging Group 
(HFAG) are also shown. 



QCDF [50] 



PQCD [54, 55] SCET [58] 



CP [92] exp 



BR{tt-K°) 

BRijT^K-) 

Acp(7r0i^-) 

BR{'K+K-) 

Acvi-i^+K-) 

BR{iT^K^) 

ylcp(vr°i?°) 



19.3 



+1.9+11.3+1.9+13.2 
-1.9- 7.8-2.1- 5.6 



n Q +0.2 +0.3 +0.1 +0.6 
'-'•^-0.3-0.3-0.1-0.5 

1 1 1 +1.8+5.8+0.9+6.9 
-1.7-4.0-1.0-3.0 

7 1 +1.7+2.0+0.8+9.0 
'■^ -1.8-2.0-0.6-9.7 



16.3 



4.5 



+2.6+9.6+1.4+11.4 
-2.3-6.5-1.4- 4.8 

+1.1+2.2+0.5+8.7 
-1.1-2.5-0.6-9.5 



7 r, +0-7 +4.7 +0.7 +5.4 
'■'-'-0.7-3.2-0.7-2.3 

q q+1. 0+1.3+0.5+3.4 

''■'^-0.8-1.6-1.0-3.3 



24.5 



+13.6 
- 8.1 



0±0 
13.91^°;° 



20.9 



+15.6 
8.3 



20.8 ± 7.9 ± 0.6 ± 0.7 24.1 ±0.7 23.1 ± 1.0 

<5 1.2 ±2.4 0.9 ±2.5 

11.3 ±4.1 ±1.0 ±0.3 12.6 ±0.5 12.8 ± 0.6 

-11±9±11±2 3.4 ±2.4 4.7 ±2.6 

20.1 ±7.4 ±1.3 ±0.6 19.6 ±0.5 19.4 ± 0.6 

-6±5±6±2 -8.9 ±1.6 -9.5 ±1.3 

9.4 ±3.6 ±0.2 ±0.3 9.5 ±0.4 10.0 ±0.6 

5±4±4±1 -9.8 ±3.7 -12 ±11 



TABLE II: Results for two-body b ^ s penguin decays to rj or r]' CP-averaged BR^s (in unit of 
10~^) and CP asymmetries (in %) in several approaches. Experimental averages from HFAG are 
also shown. 

QCDF [50] SCET I [58] SCET II [58] exp 

BR{K%') 46.5 t^^i^^ji^^g -"13^5° 63.2 ± 24.7 ± 4.2 ± 8.1 62.2 ± 23.7 ± 5.5 ± 7.2 64.9 ± 3.5 
^CP(-^V) l-Slaslai^oilol 1.1 ± 0.6 ± 1.2 ± 0.2 -2.7 ± 0.7 ± 0.8 ± 0.5 9±6 
BR{KS) l-ltlltlltoitol 2.4 ±4.4 ±0.2 ±0.3 2.3 ± 4.4 ± 0.2 ± 0.5 < 1.9 
Acp{K^r]) -9.0t|jlJ;^l|^t?| 21 ± 20 ± 4 ± 3 -18±22±6±4 
BR{K-T]') '^9Atlitllit7A^JifQ^ 69.5 ± 27.0 ± 4.3 ± 7.7 69.3 ± 26.0 ± 7.1 ± 6.3 69.711? 
Acp{K-r)') ^Ato^toitoitli -1 ± 0.6 ± 0.7 ± 0.5 0.7 ± 0.5 ± 0.2 ± 0.9 3.1 ± 2.1 
BR{K-ri) l.9toitlitli1:l-^ 2.7 ± 4.8 ± 0.4 ± 0.3 2.3 ± 4.5 ± 0.4 ± 0.3 2.2 ± 0.3 
Acp{K-ri) -18.9tl\lt\7itlit^2ii 33±30±7±3 -33 ±39 ±10 ±4 29 ±11 



complete generality of the method. 

Notice also that B — Kn data in Tab. I are perfectly reproduced in the GP approach, thus 
showing on general grounds the absence of any "/Ttt puzzle" , although specific dynamical 
assumptions may lead to discrepancies between theory and experiment [80, 97-100]. 
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TABLE III: Results for CP-averaged BR^s (in units of 10 ®) and CP asymmetries (in %) in several 
approaches for B — > PV decays. Experimental averages from HFAG are also shown. 





QCDF [50] 


PQCD [54, 55] 


CP [92] 


exp 


BR{tt-K*^) 


o f; +0.4 +1.5 +1.2 +7.7 
"-••"-0.3-1.4-1.2-2.3 


6.ot?:^ 


11.3 ±0.9 


10.7 ±0.8 




1 R +0.4 +0.6 +0.5 +2.5 
-'^•"-0.5-0.5-0.4-1.0 




-7±6 


-8.5 ±5.7 


BR^-K^K*-) 


o o +1.1 +1.0 +0.6 +4.4 
"'•''-1.0-0.9-0.6-1.4 


4.311° 


7.3 ± 0.6 


6.9 ±2.3 




o 7+2.1+5.0+2.9+41.7 
°-' -2.6-4.3-3.4-44.2 


q9+21 
'''^-28 


-2± 13 


4 ±29 


BR{tt+K*-) 


Q o +1.4 +1.3 +0.8 +6.2 
"'••'-1.2-1.2-0.8-1.6 


'-'•^-2.6 


8.5 ±0.8 


9.8 ± 1.1 


Acp{n+K*-) 


9 -1 +0.6+8.2+5.1+62.5 
"^■^ -0.7-7.9-5.8-64.2 


-60l?2 


-4± 13 


-5± 14 


BRii^^K*^) 


n 7+0.1+0.5+0.3+2.6 
"•' -0.1-0.4-0.3-0.5 




3.1 ±0.4 


o.olJi 


Acp(7r0^*0) 


10 o +4.0 +4.7+2.7+31.7 
-■-^•^ -3.2 -7.0 -4.0 -35.3 


-nil 


-11± 15 


-1±27 


BR{K^p-) 


r n +0.6 +7.0+1.5 +10.3 
"J^® -0.6 -3.3 -1.3- 3.2 


o 7+6.8 


7.8 ±1.1 


8.011:1 


AcpiK'p-) 


n Q +0.1 +0.3 +0.2 +1.6 
'-'••'-0.1-0.4-0.1-1.3 


1 ± 1 


0.02 ±0.17 


12± 17 


BR{K'p^) 


9 +0.9 +3.1 +0.8 +4.3 
^•'^-0.9-1.4-0.6-1.2 


^•l^ts 


4.15 ±0.50 


^•^"^-0.56 


Acp(K-pO) 


iQ +4.5 +6.9 +3.7 +62.7 
-'-''•'-' -5.7 -4.4 -3.1 -55.4 


71 +25 
' ^-35 


29± 10 


3lliJ 


BR{K-p+) 


7 A +1.8+7.1+1.2+10.7 
'•^-1.9-3.6-1.1- 3.5 


8.812:1 


10.2 ± 1.0 


15.311^ 


Acp{K-p+) 


o Q +1.3 +4.4 +1.9 +34.5 
•'•0-1.4 -2. 7-1. 6-32.7 


64li^ 


21 ±10 


22 ±23 


BR{K^p^) 


A +0.5 +4.0 +0.7 +6.1 
^•'^-0.5-2.1-0.7-2.1 


4.8l|i 


5.2 ±0.7 


^ 4-HO.9 
^•^-1.0 


AcpiK'p') 


7 c- +1.7+2.3+0.7+8.8 
'•^-2.1-2.0-0.4-8.7 


7+8 
'-5 


1±15 


-64 ± 46 


BR{K-u) 


q cr +1.0 +3.3 +1.4 +4.7 
''•^-1.0-1.6-0.9-1.6 


10.6l^°3^ 


6.9 ±0.5 


6.8 ±0.5 


AcpiK~uj) 


7 o +2.6 +5.9 +2.4 +39.8 
'•° -3.0 -3.6 -1.9 -38.0 


321J5 


5±6 


5±6 


BR{K^uj) 


n q +0.3 +2.8 +1.3 +4.3 
^••'-0.3-1.3-0.8-1.3 


9.8li« 


4.6 ±0.5 


5.2 ±0.7 


AcpiK'^uj) 


o 1 +2.5+3.0+1.7+11.8 
"•-"^ -2.0 -3.3 -1.4 -12.9 


q+2 
-0_4 


-5± 11 


21 ± 19 


BR{K-(P) 


A c- +0.5 +1.8 +1.9 +11.8 
^•0-0.4-1.7-2.1- 3.3 


7.8l?| 


8.39 ± 0.59 


8.30 ±0.65 


Acp{K-ct>) 


1 +0.4 +0.6 +0.5 +3.0 
^•"^ -0.5 -0.5 -0.3 -1.2 


11? 


3.0 ±4.5 


3.4 ±4.4 


BR{K^^) 


A 1 +0.4+1.7+1.8+10.6 
^•J^ -0.4-1.6-1.9- 3.0 


7 q+5.4 

' •''-1.6 


7.8 ± 0.9 


8-311:^ 


Acp(^V) 


1 7+0.4+0.6+0.5+1.4 
■■■•' -0.5-0.5-0.3-0.8 


o+l 
"3-2 


1±6 


-1±13 



We conclude that factorization approaches in general show a remarkable agreement with 
experimental data, but their predictions suffer from large uncertainties. Furthermore, QCDF 
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TABLE IV: Results for CP-averaged BR^s (in units of 10 ®) and CP asymmetries (in %) in several 
approaches for Bs PP decays. The only available experimental result is BR{Bs —>■ K~^K~) = 
(24.4 ±4.8) • IQ-^ [96]. 

QCDF [50] PQCD [95] SCET I [58] SCET 11 [58] 

BR{K+K-) 22.7l|it^|Il|°12^;} 17.0tl:lt|^l°:l 18.2 ± 6.7 ± 1.1 ± 0.5 

AMK-^K-) 4.oil:ri:°i°:ii}?:| -25.8iJiit:li?:? -6±5±6±2 

BR{K^K^) 24.7 +2.6+25.6 ig.el^;^!^"!™ 17.7 ± 6.6 ± 0.5 ± 0.6 

Acp(K°K°) o.9t°;^2lSltSilSJ < 10 

BRir^r]) 15.6l{ 14.6l|^l|n:[J 7.1 ± 6.4 ± 0.2 ± 0.8 6.4 ± 6.3 ± 0.1 ± 0.7 

^cp(w) -i.6l°ii°Jir7lli -i-6lg;ilgilgi 7.9 ±4.9 ±2.7 ±1.5 -1.1 ±5.0 ±3.9 ±1.0 
BR{-qi) 54.0 111 11^1^^:^ 39.ol?:^l?^:^lg:g 24.0 ±13.6 ±1.4 ±2.7 23.8 ±13.2 ±1.6 ±2.9 
Acp(w') o.4i|];{i[jj;[|:}i[!:^ -i.2t{5:olo:?;oi 0.04 ±0.14 ±0.39 ±0.43 2.3 ± 0.9 ± 0.8 ± 7.6 

BR{r]'ri') 4L7ltoi?7:ii1:5i?5:4 29.6l|^t^1j%°o° 44.3 ± 19.7 ±2.3 ± 17.1 49.4 ± 20.6 ± 8.4 ± 16.2 
AcpCVV) 2.1 10:510:4+0-2+i.i 2.2l0;nil0;2 0.9 ± 0.4 ± 0.6 ± 1.9 -3.7 ± 1.0 ± 1.2 ± 5.6 



and SCET cannot reproduce rate asymmetries in S ^ i^vr; this might be a hint that some 
dehcate aspects of the dynamics of penguin decays, for example rescattering and final state 
interaction phases, are not fully under control. It is then reassuring that a more general 
approach as GP can reproduce the experimental data with reasonable (but not too small) 
values of the A/mf, corrections to factorization. To quantify this statement, we report in 
Fig. 1 the results of the GP fit for Acr{B Kn) as a function of the upper bound on 
A/rrib corrections [92]. It is clear that imposing a too low upper bound, of order 10%, would 
generate a spurious tension between theory and experiment. 

For the reader's convenience, we report in Tabs. IV and V the predictions obtained in 
several approaches for BR^s and GP asymmetries of Bg penguin-dominated 6 — > s decays. 



B. Predictions for S and AS in s penguins 

Keeping in mind the results of Sec. Ill A, we now turn to the main topic of this review, 
namely our ability to test the SM using time-dependent GP asymmetries in 6 — > s penguin 
nonleptonic decays. 
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TABLE V: Results for CP-averaged BR^s (in units of 10 ^) and CP asymmetries (in %) in several 
approaches for Bg —>■ PV decays. No experimental data are available yet. 



Channel QCDF [50] PQCD [95] 

mi(T^+T^*-\ /I 1 +1.7+1.5+1.0+9.2 7/1+2.1+1.9+0.9 

nn\^J\ J\ ) "i-J^ -1.5 -1.3 -0.9 -2.3 '-^-l. 8-1. 4-0.4 

Ar.T,(K+K*-\ o 9+0.6+8.4+5.1+68.6 ^.+2.9+2.2+1.8 

j^CP\J^ ) ^•^-0.7-8.0-5.9-71.0 ^"J-0-2.4-3.0-1.3 

FtJilK^ Q +0-4 +1-5 +1-3 +10.4 q i +3.2+2.6+0.0 

nnyix J\ ) ^-^-0.4-1.4-1.4- 2.8 ^■-'^-2.2-1.5-0.0 

Rn(T<r-K*+\ c; +1.3 +5.0 +0.8 +14.2 ^ c-+l. 2+3.3+0.0 

Dn\^J\ J\ ) O-O_i,4_2.6-0.7- 3.6 0-^-1.2-1.8-0.1 

Ar.-r.(K-K*+\ _o 1 +1.0+3.8+1.6+47.5 f^o 0+5.2+8.0+5.1 

^CPl-f^ J\ ; <J-J- -1.1 -2.6 -1.3 -45.0 ""^-^-5.8-10.2-2.6 

RJi(K^K*^\ 4 9+0.4+4.6+1.1+13.2 c Q+0-9+2.8+0.0 

nnyjx J\ ) "±-^_o.4-2.2-0.9- 3.2 ^-^-1.1-1.8-0.0 

R j^f^, A n ni 9 +0.005 +0.010 +0.028 +0.025 ^ 1 n+0-02+o.03+o.oo 



AcpM 3.2 



^-0.02-0.01-0.00 

+6.1+15.2+0.3 
-3.9-11.2-0.1 



UTK^', ,\ n no/l +0-011 +0-028 +0.077 +0.042 n. ^^.+0. 23+0. 22+0. 01 

orLyijLO) U-*-*^^ -0.009 -0.006 -0.010 -0.015 '-'-"'J-0.18-0.21-0.03 

Rnl'nM 19+0.02+0.95+0.54+0.32 i O+0.5+0.1+0.0 

JDJX\ll(p) "--'-^-0.02-0.14-0.12-0.13 ^-°-0.5-0.2-0.0 

(r.AA 0^+2.0+30.1+14.6+36.3 r, 1 +0-2+2.3+0.0 

^CPlWj -f5-4_2.i -71.2-44.7 -59.7 -^-^-0.4-1.4-0.0 

Rpr^iU^ n nc^ +0-01 +1.10 +0.18 +0.40 0^+1.2+0.4+0.0 

nnyrjcp) "-"J-J -O.Ol -0.17-0.08 -0.04 •J-'-'-o. 9-0. 4-0.0 

Ar.r.t'n'rh') fio 9+15.9+132.3+80.8+122.4 i 9+O.I+O.4+O.I 

^CP{y 9) ""J^-^ -10.2- 84.2-46.8- 49.9 -'--^-0.0-0.6-0.1 



Starting from Eq. (5), we write down the expression for 5*^ as follows: 



sin(2(/3, + (f)M)) + 1 


\i^f\ 


2 sin(2((/)M + 7)) + 2 Re sm((3s + '^(f^u + 7) 


1 + 




|2 + 2Rerircos(/?, -7) 



(7) 



where rp = |KsK&|/|V^.^t6| x E(71 + P^™)/Y.Pi with = for pure penguin channels. 
Since the angle /3s is small and very well known (/3s — (2.1 ± 0.1)°), the problem is then 
reduced to the evaluation oi Kp = ^{Ti + P[^^^'^) / J2Pi for each channel (notice that Ti = 
for pure penguin channels). Factorization methods have been used to provide estimates of 
Kp, Sp and ASp for 6 ^ s channels. The latter are reported in Tables VI and VII. A 
few remarks are important. First of all, the evaluation of p9^^ rehes on the factorization 
of penguin contractions of charm and up quarks, which is debatable even in the infinite 
mass limit. In addition to that, in factorization P^™^ has a perturbative loop suppression 
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Uppi^r Value 0(A^^m^ Upper Value OlA^-^m^ Upper Value 0(A^^m^ 



FIG. 1: j4cp values for B Kn in the GP approach [92], obtained varying 0{Aqcd /f^b) con- 
tributions in the range [0, UV], with the upper value UV scanned between zero and one (in units 
of the factorized emission amplitude). For comparison, the experimental 68% (95%) probability 
range is given by the dark (light) band. 

TABLE VI: Predictions for S parameters in % for B decays. Experimental averages from HFAG 
are also shown. 



PQCD [54, 55] SCET I [58] SCET II [58] GP [92] exp 





74: 


t-2 


80 ± 2 ± 2 ± 1 


74.3 ±4.4 33 ±21 


S-q'Ks 






70.6 ± 0.5 ± 0.6 ± 0.3 71.5 ± 0.5 ± 0.8 ± 0.2 


70.9 ±3.9 61 ±7 


SrjKs 






69 ± 15 ± 5 ± 1 79 ± 14 ± 4 ± 1 




S^Ks 


71: 


-1 




71.5 ±8.7 39 ±18 


SpOKs 


50l 


10 
6 




64 ±11 20 ± 57 


SujKs 


84: 


^3 
~7 




75.7 ± 10.3 48± 24 



so that it is likely to be dominated by power corrections. Furthermore, the contribution of 
Tj and P^™ is particularly difficult to estimate for r] and 77' channels. Last but not least, 
the determination of the sign of ^Sp relies heavily on the determination of the sign of Re 
Kp. If P^™ is dominated by power corrections, there is no guarantee that the sign given by 
the perturbative calculation is correct. 

With the above caveat in mind, from Tables VI and VII we learn that: 

• Experimentally there is a systematic trend for negative AS. This might be a hint of 
the presence of new sources of CP violation in the b s penguin amplitude. 

• The experimental uncertainty is dominant in all channels. In addition to that, the GP 
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TABLE VII: Predictions for AS parameters in % for B decays. Experimental averages from HFAG 
are also shown. 





QCDF [101] 


SCET I [58] 


SCET II [58] 


GP [92] 


exp 


AS^rOKs 


7+5 
'-4 


7.7 ± 2.2 ± 1.8 ± 1 




2.4 ±5.9 


-35 ±21 






-1.9 ±0.5 ±0.6 ±0.3 


-1.0 ±0.5 ±0.8 ±0.2 


-0.7 ±5.4 


-7 ±7 


ASrjKs 


lOl^i 


-3.4 ± 15.5 ±5.4 ± 1.4 


7.0 ± 13.6 ±4.2 ± 1.1 






ASrt>Ks 








0.4 ±9.2 


-29 ± 18 


ASpOKs 


o+ 8 






-6.2 ±8.4 


-48 ± 57 


AS^^Ks 


1311 






5.6 ± 10.7 


-20 ± 24 



estimate of the theoretical uncertainty, which is certainly conservative, can be reduced 
with experimental improvements on BRs and CP asymmetries. 

• As discussed in Sec. II, the theoretical uncertainty estimated from first principles is 
much smaller for pure penguin decays such as 5 — > 4>Ks than for penguin-dominated 
channels. 

• In the mo del- independent GP approach, the theoretical uncertainty is smaller for 
B — ^ TT^Kg because the number of observables in the B — >■ Ktt system is sufficient to 
constrain efficiently the hadronic parameters. This means that the theoretical error 
can be kept under control by improving the experimental data in these channels. 
On the other hand, the information on S — > (pKg is not sufficient to bound the 
subleading terms and this results in a relatively large theoretical uncertainty that 
cannot be decreased without additional input on hadronic parameters. Furthermore, 
using SU (3) to constrain AS^^Ks difficult because the number of amplitudes involved 
is very large [64, 85-87]. 

The ideal situation would be represented by a pure penguin decay for which the infor- 
mation on p9™ is available with minimal theoretical input. Such situation is realized by 
the pure penguin decays Bg i^o(*)^o(*). An upper bound for the Pj^^^ entering this 
amplitude can be obtained from the 5'C/(3)-related channels B^ — > K^(*') K°(*') . Then, even 
adding a generous 100% SU{3) breaking and an arbitrary strong phase, it is possible to have 
full control over the theoretical error in AS [94]. 
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TABLE VIII: Predictions for S parameters for Bs decays. 



PQCD [95] SCET I [58] SCET II [58] 
-0.16 ±0.41 ±0.33 ±0.17 

0.82 ± 0.32 ± 0.11 ± 0.04 0.63 ± 0.61 ± 0.16 ± 0.08 

0.38 ± 0.08 ± 0.10 ± 0.04 0.24 ± 0.09 ± 0.15 ± 0.05 
0.19 ±0.04 ±0.04 ±0.01 

0.45 ± 0.14 ± 0.42 ± 0.30 0.38 ± 0.20 ± 0.42 ± 0.37 

-0.026 ± 0.040 ± 0.030 ± 0.014 -0.077 ± 0.061 ± 0.022 ± 0.026 

0.041 ± 0.004 ± 0.002 ± 0.051 0.015 ± 0.010 ± 0.008 ± 0.069 

0.049 ± 0.005 ± 0.005 ± 0.031 0.051 ± 0.009 ± 0.017 ± 0.039 





KsTT^ 


n /ic+0. 14+0. 19+0.02 
^•^"-'-0.13-0.20-0.04 




KsV 


p. oi +0.05+0.16+0.02 
"••Jj--0.05-0.17-0.03 




Ksi 


r, 79+0.02+0.04+0.00 
•-"^ '^-0.02-0.03-0.00 


B's^ 


K-K+ 


n oQ+0.04+0.04+0.02 
^•^°-0.04-0.03-0.01 




Tr% 


n nn+0-03+0.09+0.00 

^•"^-0.02-0.10-0.01 




rjrj 


n n'^+o-oo+o.oi+o.oo 
"•"•J-o.oo-o.oi-o.oo 




rji]' 


n r,/i +0.00+0.00+0.00 

'^•"^-0.00-0.00-0.00 




r]'r]' 


n n/i +0.00+0.00+0.00 
"•U^-0.00-0.00-0.00 




corj 


n 07+0.00+0.04+0.00 
"•"'-0.01-0.11-0.00 




corj' 


n 1 Q+0.01+0.04+0.01 

"•^^-0.01-0.04-0.03 




4>v 


rv -in+0-01+0.04+0.01 

"•-■-"-0.01-0.03-0.00 




4>v' 


n nn+0.00+0.02+0.00 
"•""-0.00-0.02-0.00 




Ks<t> 


-0.72 



For the reader's convenience, we report in Tab. VIII the predictions for the S coefficient 
of the time-dependent CP asymmetry for several Bg penguin-dominated decays. 

Before closing this Section, let us mention non-resonant three-body B decays such as 
B KsTi^Ti'^, B KsKgKs or B K^K^Kg. In this case, a theoretical estimate of Kp 
is extremely challenging, and using SU{3) to constrain kf is difficult because of the large 
number of channels involved [85]. Nevertheless, they are certainly helpful in completing the 
picture of CP violation in 6 — > s penguins. 

To summarize the status of 6 ^ s penguins in the SM, we can say that additional 
experimental data will allow us to establish whether the trend of negative AS shown by 
present data really signals the presence of NP in 6 — > s penguins. Theoretical errors are not 
an issue in this respect, because the estimates based on factorization can in most cases be 
checked using the GP approach based purely on experimental data. Bg decays will provide 
additional useful channels and will help considerably in assessing the presence of NP in 
6 — > s penguins. 
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IV. CP VIOLATION IN 6 ^ s PENGUINS BEYOND THE SM 

We have seen that there is a hint of NP in CP-violating 6 — > s hadronic penguins. In 
this Section, we would like to answer two basic questions that arise when considering NP 
contributions to these decays: 

1. What are the constraints from other processes on new sources of CP violation in 6 — > s 
transitions? 

2. Are NP contributions to b ^ s transitions well motivated from the theoretical point 
of view? 

We consider here only model-independent aspects of these two questions, and postpone 
model-dependent analyses to Section V. 



A. Model-independent constraints on s transitions 

The last year has witnessed enormous progress in the experimental study oi b ^ s tran- 
sitions. In particular, the TeVatron experiments have provided us with the first information 
on the Bs — Bs mixing amplitude [102], which can be translated into constraints on the 
AB — AS — 2 effective Hamiltonian. In any given model, as we shall see for example in 
Sec. V, these constraints can be combined with the ones from 6 — > 57 and b — > sl'^l~ decays 
to provide strong bounds on NP effects in 6 — > s hadronic penguins. 

Let us now summarize the presently available bounds on the Bg — Bg mixing amplitude, 
following the discussion of rcf. [103]. General NP contributions to the Ai? = AiS = 2 effective 
Hamiltonian can be incorporated in the analysis in a model-independent way, parametrizing 
the shift induced in the mixing frequency and phase with two parameters, Cb^ and 0^^, 
equal to 1 and in the SM [104-108]: 

As for the absorptive part of the Bg — Bg mixing amplitude, which is derived from the 
double insertion of the AB — 1 effective Hamiltonian, it can be affected by non-negligible 
NP effects in AB — 1 transitions through penguin contributions. Following refs. [4, 5], we 
thus introduce two additional parameters, Cf^" and 0^^", which encode NP contributions 
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FIG. 2: Constraints on the (pB^ vs. Cb^ plane [103]. Darker (lighter) regions correspond to 68% 
(95%) probability. 

to the penguin part of the AB = 1 Hamiltonian in analogy to what Cb^ and (pBs do for the 
mixing amphtude. 

The available experimental information is the following: the measurement of Arris [102], 
the semileptonic asymmetry in Bg decays A|l and the dimuon asymmetry Aqu from D0 [109, 
110], the measurement of the Bg lifetime from flavor-specific final states [111-116], the 
determination of AFs/F^ from the time-integrated angular analysis of Bs J/'4"P decays 
by CDF [117], the three-dimensional constraint on F^, AF<j, and Bg-Bg mixing phase (pg 
from the time-dependent angular analysis of Bg J/'ipfj) decays by D0 [118]. 

Making use of this experimental information it is possible to constrain Cbs and (pSs [5, 
103, 119-122]. The fourfold ambiguity for (pSs inherent in the untagged analysis of ref. [118] 
is somewhat reduced by the measurements of A^^ and Asl [123], which prefer negative values 
of (pBs- The results for Cb^ and 0^^, obtained from the general analysis allowing for NP in 
all sectors, are [103] 

= 1.03 ±0.29, </)B^ = (-75±14)°U (-19±11)°U(9±10)°U(102±16)° . (9) 

Thus, the deviation from zero in is below the la level, although clearly there is still 
ample room for values of (pBs very far from zero. The corresponding p.d.f. in the Cb,-4'Bs 
plane is shown in Fig. 2. 
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The experimental information on 6 — > 57 and b — > si'^i" decays [124-130] can also 
be combined in a model-independent way along the lines of refs. [131-134]. In this way, 
it is possible to constrain the coefficients of the 6 — >^ S7, 6 — > 57* and b ^ sZ vertices, 
which also contribute to 6 ^ s hadronic penguins. It turns out that order-of-magnitude 
enhancements of these vertices are excluded, so that they are unlikely to give large effects 
in 6 — > s nonleptonic decays. On the other hand, the b ^ sg vertex is only very weakly 
constrained, so that it can still give large contributions to 6 — > s hadronic penguins. Finally, 
the information contained in Eq. (9) can be used to constrain NP effects in & — > s hadronic 
decays only within a given model, since a connection between AS — 2 and AS = 1 effective 
Hamiltonians is possible only once the model is specified. We shall return to this point in 
Sec. V. 



B. Theoretical motivations for NP in 6 — ^ s transitions 

We now turn to the second question formulated at the beginning of this Section, namely 
whether on general grounds it is natural to expect NP to show up in 6 — > s transitions. 
The general picture emerging from the generalized Unitarity Triangle analysis performed in 
ref. [4, 5, 103] and from the very recent data on D — D mixing [135-138] is that no new 
sources of CP violation arc present in Bd-, K and D mixing amplitudes. Conversely, large NP 
contributions to s ^ dg, b ^ dg and b ^ sg transitions are not at all excluded. Therefore, 
although the idea of minimal flavor violation is phenomenologically appealing [21, 139-144], 
an equally possible alternative is that NP is contributing more to AF — 1 transitions than 
to AF = 2 ones. Within the class of AF = 1 transitions, (chromo)-magnetic vertices are 
peculiar since they require a chirality flip to take place, which leads to a down-type quark 
mass suppression within the SM. On the other hand, NP models can weaken this suppression 
if they contain additional heavy fermions and/or additional sources of chiral mixing. In this 
case, they can lead to spectacular enhancements for the coefficients of (chromo)-magnetic 
operators. Furthermore, if the relevant new particles are colored, they can naturally give a 
strong enhancement of chromomagnetic operators while magnetic operators might be only 
marginally modifled [145]. The electric dipole moment of the neutron puts strong constraints 
on new sources of CP violation in chirality-flipping flavor-conserving operators involving light 
quarks, but this does not necessarily imply the suppression of flavor-violating operators. 
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especially those involving b quarks. Therefore, assuming that NP is sizable in hadronic 
6 — > s penguins is perfectly legitimate given the present information available on flavor 
physics. 

From a theoretical point of view, a crucial observation is the strong breaking of the 
SM SU{3)^ flavor symmetry by the top quark Yukawa coupling. This breaking necessarily 
propagates in the NP sector, so that in general it is very difficult to suppress NP contributions 
to CP violation in b decays, and these NP contributions could be naturally larger in 6 — > s 
transitions than in & — > d ones. This is indeed the case in several flavor models (see for 
example Ref. [146]). 

Another interesting argument is the connection between quark and Icpton flavor violation 
in grand unified models [147-150]. The idea is very simple: the large flavor mixing present 
in the neutrino sector, if mainly generated by Yukawa couplings, should be shared by right- 
handed down-type quarks that sit in the same SU{5) multiplet with left-handed leptons. 
Once again, one expects in this case large NP contributions to 6 — > s transitions. 

We conclude that the possibihty of large NP effects in 6 — > s penguin hadronic decays 
is theoretically well motivated on general grounds. The arguments sketched above can of 
course be put on firmer grounds in the context of specific models, and we refer the reader 
to the rich literature on this subject. 

V. SUSY MODELS 

Let us now focus on SUSY and discuss the phenomenological effects of the new sources 
of flavor and CP violation in 6 ^ s processes that arise in the squark sector [151-173]. 
In general, in the MSSM squark masses are neither flavor-universal, nor are they aligned 
to quark masses, so that they are not flavor diagonal in the super-CKM basis, in which 
quark masses are diagonal and all neutral current vertices are flavor diagonal. The ratios of 
off-diagonal squark mass terms to the average squark mass define four new sources of flavor 
violation in the 6 — > s sector: the mass insertions {62^) ab, with A,B — L,R referring to 
the hehcity of the corresponding quarks. These 5's are in general complex, so that they also 
violate CP. One can think of them as additional CKM-type mixings arising from the SUSY 
sector. Assuming that the dominant SUSY contribution comes from the strong interaction 
sector, i.e. from gluino exchange, all FCNC processes can be computed in terms of the 
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SM parameters plus the four (5's plus the relevant SUSY parameters: the gluino mass m^, 
the average squark mass rriq, tan/3 and the /i parameter. The impact of additional SUSY 
contributions such as chargino exchange has been discussed in detail in Ref. [166]. We 
consider only the case of small or moderate tan f3, since for large tan f3 the constraints from 
Bs —>■ and Arris preclude the possibility of having large effects in 6 — > s hadronic 

penguin decays [163, 174-179]. 

Barring accidental cancellations, one can consider one single 5 parameter, fix the SUSY 
masses and study the phenomenology. The constraints on (5's come at present from B — > 
Xg^, B ^ Xsl'^l~ and from the Bg — Bg mixing amplitude as given in Eq. (9). We refer the 
reader to rcfs. [180-182] for all the details of this analysis. 

Fixing as an example rrig = rriq = = 350 GeV and tan/? = 3 or 10, one obtains 
the constraints on 5's reported in Figs. 3-5 [181, 182]. We plot in light green the allowed 
region considering only the constraint from the Cb^ vs. p.d.f. of Fig. 2, in hght blue 
the allowed region considering only the constraint from b — > s£'^£~ and in violet the allowed 
region considering only the constraint from b ^ s^. The dark blue region is the one selected 
imposing all constraints simultaneously. 

Several comments are in order at this point: 

• Only (523)ll,lr generate amphtudes that interfere with the SM in rare decays. There- 
fore, the constraints from rare decays for (523)rl,rr are symmetric around zero, while 
the interference with the SM produces the circular shape of the B — > Xg'y constraint 
on (523)ll,lr- 

• We recall that LR and RL mass insertions generate much larger contributions to the 

(chro mo) magnetic operators, since the necessary chirality flip can be performed on the 
gluino line (oc rrig) rather than on the quark line (oc m^). Therefore, the constraints 
from rare decays are much more effective on these insertions, so that the bound from 
Bg — Bg has no impact in this case. 

• The //.tan/5 flavor-conserving LR squark mass term generates, together with a flavor 
changing LL mass insertion, an effective (523)lr that contributes to B ^ Xgj. For 
positive (negative) /x, we have ((^23)lr +(— )(<^23)ll and therefore the circle deter- 
mined by S — > Xg^ in the LL and LR cases lies on the same side (on opposite sides) 
of the origin (see Figs. 3 and 4). 
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FIG. 3: Allowed region in the Ref ) -Imf ) plane. In the plots on the left (right), negative 

v / LL \ / LL 

(positive) fi is considered. Plots in the upper (lower) row correspond to tan/? = 3 (tan/? = 10). 
See the text for details. 



• For tan/3 = 3, we see from the upper row of Fig. 3 that the bound on ((523)ll from 
Bg — Bg mixing is competitive with the one from rare decays, while for tan /5 = 10 rare 
decays give the strongest constraints (lower row of Fig. 3). The bounds on all other 
5's do not depend on the sign of and on the value of tan (3 for this choice of SUSY 
parameters. 

• For LL and LR cases, B —>■ Xg'-f and B Xgl^l^ produce bounds with different 
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shapes on the Re 5 - Im 5 plane (violet and light blue regions in Figs. 3 and 4), so 
that applying them simultaneously a much smaller region around the origin survives 
(dark blue regions in Figs. 3 and 4). This shows the key role played by rare decays in 
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FIG. 6: Allowed region in the Ke(62^] -ImfJis) plane. In the plots on the left (right), 

negative (positive) /i is considered. Plots in the upper (lower) row correspond to tan/3 = 3 (tan/3 = 



10). See the text for details. 



constraining new sources of flavor and CP violation in the squark sector. 

• For the RR case, the constraints from rare decays are very weak, so that the only 
signiflcant bound comes from — Bg mixing. 

• If (5f3)LL and (523)rr insertions are simultaneously nonzero, they generate chirality- 
breaking contributions that are strongly enhanced over chirality-conserving ones, so 
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that the product (<523)ll(^23)rr severely bounded. In Fig. 6 we report the allowed 
region obtained in the case (523)ll = ('^23)rr- ^'^^ (^23)ll 7^ ('^23)rR' ^^is constraint can 
be interpreted as a bound on (523)ll(<^23)rr- We observe a very interesting interplay 
between the constraints from rare decays and the one from Bg — Bg mixing. Increasing 
tan /3 from 3 to 10, the bound from rare decays becomes tighter, but Bg — Bg mixing 
still plays a relevant role. 

• All constraints scale approximately linearly with squark and gluino masses. 




FIG. 7: Probability density functions for S^Ks^ ^-k^Ks^ '^v'Ks ^^^^ Si^k^ induced by (523)ll- 
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FIG. 8: Probability density functions for 5^/^^, St^ox^, S^iKs ^ujKs induced by ((523)lr- 

Having determined the p.d.f 's for the four 5's, we now turn to the evaluation of the time- 
dependent CP asymmetries. As we discussed in Sec. II, the uncertainty in the calculation of 
SUSY effects is even larger than the SM one. Furthermore, we cannot use the GP approach 
since to estimate the SUSY contribution we need to evaluate the hadronic matrix elements 
explicitly. Following ref. [180], we use QCDF, enlarging the range for power-suppressed 
contributions to annihilation chosen in Ref. [50] as suggested in Ref. [42]. We warn the 
reader about the large theoretical uncertainties that affect this evaluation. 

In Figs. 7-10 we present the results for S^k^i St,ok^, Sr,'Ks and Sujk^- They do not show a 
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FIG. 9: Probability density functions for 5^/^^, St^ox^, S^iKs ^ujKs induced by ((523)rl- 



sizable dependence on the sign of or on tan (3 for the chosen range of SUSY parameters. 
We see that: 

• deviations from the SM expectations are possible in all channels, and the present 
experimental central values can be reproduced; 

• they are more easily generated by LR and RL insertions, due to the enhancement 
mechanism discussed above. 

• As noticed in refs. [183, 184], the correlation between ASpp and ASpv depends on 
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FIG. 10: Probability density functions for S^Ks-, S-j^o^^, Sy'Ks and S^jK^ induced by ('^23)rr- 



the chirality of the NP contributions. For example, we show in Fig. 11 the correlation 
between ASKs<t> and ASks-kO for the four possible choices for mass insertions. We 
see that the ASkscP and ASk^ttO are correlated for LL and LR mass insertions, and 
anticorrelated for RL and RR mass insertions. 

An interesting issue is the scaling of SUSY effects in AS* with squark and gluino masses. 
We have noticed above that the constraints from other processes scale linearly with the 
SUSY masses. Now, it turns out that also the dominant SUSY contribution to AS", the 
chromomagnetic one, scales linearly with SUSY masses as long as rrig ~ rriq ~ /i. This 
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FIG. 11: Correlation between S^Ks S^^ox^ for LL, LR, RL and RR mass insertions. 



means that there is no decouphng of SUSY contributions to AS* as long as the constraint 
from other processes can be saturated for 6 < 1. From Figs. 3-5 we see that the bounds on 
LL and RR mass insertions quickly reach the physical boundary at 5 = 1, while LR and RL 
are safely below that bound. Chirality flipping LR and RL insertions cannot become too 
large in order to ensure the absence of charge and color breaking minima and unbounded 
from below directions in the scalar potential [185]. However, it is easy to check that the 
flavor bounds given above are stronger for SUSY masses up to (and above) the TeV scale. 
We conclude that LR and RL mass insertions can give observable effects for SUSY masses 
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within the reach of LHC and even above. This is shown exphcitly in Figs. 12 and 13, where 
we present the p.d.f. for S^Ks^ 3^,0^^, Srj'Ks and S^^Ks ^ot SUSY masses of 1 TeV. 







FIG. 12: Probability density functions for S^pXs^ S^ox^, S^j/Rs ^iid S^^Ks induced by (^23)lr for 
rrig = rriq = = 1 TeV. 



VI. NON-SUSY MODELS 

In general, one expects sizable values of AS* in all models in which new sources of CP 
violation are present in 6 — s penguins. In particular, models with a fourth generation, both 
vectorlike and sequential, models with warped extra dimensions in which the flavor structure 
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of the SM is obtained using localization of fermion wave functions, and models with Z' gauge 
bosons can all give potentially large contributions to 6 — s penguins [186-190]. 

In any given NP model, it is possible to perform a detailed analysis along the lines of 
Sec. V, considering the constraints from — mixing and from rare B decays, plus the 
constraints from all other sectors if they are correlated with b s transitions in the model. 
On general grounds, the dominant contributions to b —>■ s hadronic decays are expected to 
come from electroweak or chromomagnetic penguins. The correlation between the induced 
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ASpp and ASpv can give a handle on the chirahty of the NP-generated operators. NP effects 
in electroweak penguin contributions are in general correlated with effects in 6 — > si'^i", 
in b ^ S7 and possibly in Z ^ bb. Depending on the flavor structure of NP, other effects 
might be seen in K ^ ttuu or in e'/e. NP effects in the chromomagnetic penguin might also 
show up in b s'f, in Bg — Bg mixing and, if there is a correlation between the B and K 
sectors, in e'/e. 

VII. CONCLUSIONS AND OUTLOOK 

We have reviewed the theoretical status of hadronic b ^ s penguin decays. We have 
shown that, in spite of the theoretical difficulties in the evaluation of hadronic matrix el- 
ements, in the SM it is possible to obtain sound theoretical predictions for the coefficient 
Sp of time- dependent CP asymmetries, using either models of hadronic dynamics or data- 
driven approaches. Experimental data show an interesting trend of deviations from the SM 
predictions that definitely deserves further theoretical and experimental investigation. 

Prom the point of view of NP, the recent improvements in the experimental study of other 
b ^ s processes such as Bg — Bg transitions or b ^ and b s£'^i~ have considerably 
restricted the NP parameter space. However, there are still several NP models, in particular 
SUSY with new sources of 6 s mixing in squark mass matrices, that can produce deviations 
from the SM in the ballpark of experimental values. In any given model, the study of 
hadronic 6 — > s penguins and of the correlation with other FCNC processes in B and K 
physics is a very powerful tool to unravel the fiavor structure of NP. 

Any NP model with new sources of CP violation and new particles within the mass reach 
of the LHC can potentially produce sizable deviations from the SM in b ^ s penguins. It 
will be exciting to combine the direct information from the LHC and the indirect one from 
fiavor physics to identify the physics beyond the SM that has been hiding behind the corner 
for the last decades. In this respect, future facilities for B physics will provide us with an 
invaluable tool to study the origin of fermion masses and of fiavor symmetry breaking, two 
aspects of elementary particle physics that remain obscure in spite of the theoretical and 
experimental efforts in ffavor physics. 
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